Marginal Fermi Liquid in a lattice of three-body bound-states. 
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We study a lattice model for Marginal Fermi liquid behavior, involving a gas of electrons coupled 
to a dense lattice of three-body bound-states. The presence of the bound-states changes the phase 
space for electron-electron scattering and induces a marginal self-energy amongst the electron gas. 
When the three-body bound-states are weakly coupled to the electron gas, there is a substantial 
window for marginal Fermi liquid behavior and in this regime, the model displays the presence 
of two relaxation times, one linear, one quadratic in the temperature. At low temperatures the 
bound-states develop coherence leading to a cross-over to conventional Fermi liquid behavior. At 
strong-coupling, marginal Fermi liquid behavior is pre-empted by a pairing or magnetic instability, 
and it is not possible to produce a linear scattering rate comparable with the temperature. We 
discuss the low temperature instabilities of this model and compare it to the Hubbard model at 
half-filling. 
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I. INTRODUCTION 

The concept of a marginal Fermi liquid (MFL) was in- 
vented as a phenomenological description of the asymp- 
totic properties of ofJiigh temperature superconductors 
in their normal stateEl. These include a linear electrical 
resistivity p ~ T over at least 2 decades in tempera- 
ture and a quasi-particle scattering rate which may be 
proportional to frequency up to energies as high as 750 
meV.Bu Even in underdoped cuprate superconductors, 
optical data indicates that marginal behavior develops 
at scales above the spin gap.Q We shall take a marginal 
Fermi liquid to be a system of fermions with an inelastic 
scattering rate of the form 



r oc max[| lo |, T] 



(1) 



where there is no significant momentum dependence of T. 
In the cuprates, the constant of proportionality is of order 
unity. Analyticity then ensures that the appropriate self- 
energy takes the form 



~ cjlnmax[| to \,T]. 



(2) 



Many proposals have been made to explain the origin of 
this unusual behavior. One prevalent idea, is that the 
marginal Fermi liquid behavior derives from scattering 
off a soft bosanic mode. This ideajmderpins the Van 
Hove scenariou, gauge theory modelaj and the quantum- 
critical scattering^ picture of the cuprates. By contrast, 
AndersonB proposes that the cuprate metal is a fully 
developed Luttinger liquid, with power-law self-energies 
which have been mis-identified as a logarithm. Since the 
soft-mode theories furnish an electron self-energy which is 
strongly momentum dependent, none of these proposals 
actually gives rise to a marginal Fermi liquid as originally 
envisaged. 



In this paper, we return to the original proposal, ask- 
ing whether a marginal Fermi liquid can form in a dense 
electronic system. We pursue an early speculation, due to 
Ruckenstein and VarmaQ, that the marginal self-energy 
might derive from the scattering of conduction electrons 
off a dispersionless localized bound-state $ at the Fermi 
energy, giving rise to the interaction: 



Hi, 



A(c+c iT c+$. 



h.c. 



(3) 



The presence of such states precisely at the Fermi surface 
would mean that the three-particle phase space grows 
linearly with energy. Inelastic scattering into the local- 
ized bound-state then leads to a marginal self-energy in 
leading order perturbation theory, as illustrated in Fig. 
1. (The localised state is represented by the dash line.) 
Since the hypothetical object at the Fermi surface scat- 
ters electrons in triplets, Ruckenstein and Varma identi- 
fied it as a three-body bound state. 




FIG. 1. Marginal self-energy diagram for the band elec- 
trons. Notation: thick lines denote band electron propaga- 
tors, dashed lines denote the localised state propagator. 



A great difficulty with this picture is that it cannot be 
made self-consistent. At the same level of perturbation 
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$(i) (o)^( 2 ) (o)$( 3 ) (0) = A$(0), 



(6) 



where A is the amplitude for forming this pole $. The 
residual interaction with the bulk spin degrees of freedom 
in the low energy world then gives rise to a vertex of 
the form A^ 1 )^ 2 ^ 3 ^. The challenge here is to see if 
such a mechanism could be generalised to a more realistic 
lattice model. The work described below is an attempt 
to make a first step in this direction. 



FIG. 2. Self-energy for the localised state. 



II. CONSTRUCTION OF MODEL 



theory that furnishes a marginal self-energy, the three- 
body state must scatter off the electrons to produce a 
self-energy term of the form shown in Fig 2. This self- 
energy correction inevitably moves the resonance of the 
three-body bound-state away from the Fermi energy, in- 
troducing an unwanted scale into the problem and caus- 
ing the singular scattering to disappear. 

An unexpected resolution of this problem recently ap- 
peared in the .context of the single impurity two-channel 
Kondo modelll2l. Marginal Fermi liquid behavior docs 
develop in this model, and the mechanism by which it 
occurs is remarkably close to the original three-body 
bound-state proposal, with one critical difference: in the 
two-channel Kondo model, the three-body bound-state 
formed at the impurity site carries no internal quantum 
numbers (spin and charge). The associated bound-state 
fcrmion is represented by a Hcrmitian operator: 

$ = (4) 

This type of fermion is known as a Majorana fcrmion0 0. 
The effective action for such a field must have the form: 



S = /"e^$(-u;)(w - E(o/))$(w). 



(5) 



Since $(— lu)&(lj) = — $(cj)$(— lu) due to the Grassma- 
nian nature of $, T,(uj) will be an odd function of fre- 
quency, so that S(0) = 0. In other words, the particle- 
hole symmetry of the $ field guarantees that its energy 
is pinned to the Fermi surface. 

In the two-channel Kondo model (after dropping 
the charge_degrees of freedom, thanks to spin-charge 
decouplingta) , the total spin S of the two conduction 
channels at the impurity site can be written in the form: 
S = -|#(0) x *(0) where * = (^(1)^(2)^(3)) is a 
triplet of Majorana fermions. These three high energy 
degrees of freedom bind at the impurity site to form the 
localised three-body bound-state $(0), represented as the 
contraction of the three fermions: 



*Such an object can always be constructed as a linear com- 
bination of two charged fermions $ = + a )• 



We now use these ideas to motivate a simple lattice 
model of a marginal Fermi Liquid. First note that the 
Hubbard model at half-fillirutcan be rewritten in a Majo- 
rana fermion representation^, by the following two steps. 
The Hubbard model is: 



H 



Hubbard 



+ C/^(cJ TCiT -l/2)(c] iCii -l/2). (7) 

i 

In the first step, we assume that the lattice has a bipartite 
structure, and do a gauge transformation on electron op- 
erators belonging to one sublattice of the bipartite lattice: 
c,; CT — > —ici„. Also let the hopping term connect a sub- 
lattice A site only to a sublattice B site, and not to other 
sublattice A sites. Then the interaction is unchanged, 
but the kinetic energy becomes: 



H k .e. =itJ2 (A+a,* *,* + Hm -)- 



(8) 



Next, rewrite the electron operators using Cj = ("J^ 1 -* — 
i^),Cii = -^(# (3) + i* (0) ) where the *^'s are Ma- 
jorana fermions, and we get: 



(9) 



(Note the sign change for the interaction.) In Majo rana 



representation, we see explicitly the 50(4) symmetryliiE 2 . 
of the Hubbard model at half-filling. 

The crucial generalisation of this paper is to break the 
50(4) symmetry down to 0(3), setting: 



t^t„ 



t 

to 



a = 1,2,3 
a = 0. 



(10) 



With = ^\ 



(0) 



i,c i,c a—1 

-c/V$ (0) * (1) * (2) * (3) . (li) 
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When tg = 0, this Hamiltonian describes a lattice of lo- 
calised three-body bound-states $i coupled to the con- 
tinuum. 

This toy model provides a simple system to study the 
properties of these pre-formed bound-states at each site. 
A microscopic model would provide an explanation of 
the origin of this symmetry breaking field leading to for- 
mation of these bound-states. In this work, we do not 
address this issue, but instead, ask whether the single 
impurity marginal Fermi Liquid mechanism survives in 
the lattice, in this model of reduced symmetry. 

We shall show that the main physical effect of the lat- 
tice (in the absence of a broken symmetry phase) is that 
the previously localised mode $i can now move from 
site to site via virtual fluctuations into the Fermi Sea: 
see Fig. 3, thereby providing the lattice coherence energy 
scale to ~ U 2 /t, below which this marginal Fermi Liq- 
uid reverts to a Fermi Liquid. Unfortunately, this means 
that the marginal Fermi liquid phenomenology does not 
persist to large U, which makes the application of our 
model to the cuprates rather problematic. 




FIG. 3. Leading order diagram that generates a dispersion 
for the fermion. Note that there are no arrows on the 
propagator lines, as all fermions are represented by Hermitian 
operators. 

To probe the marginal fermi liquid phenomenology, 
we have calculated the "optical conductivity". This 
cannot be the ordinary electrical current, because the 
Hamiltonian of Eqn.15 does not conserve particle num- 
ber unless to — t (this is most easily seen with orig- 
inal electron operators c, c + ). The model does have 
the 0(3) symmetry, so there is the conserved quantity 

s a - -sEjEUi*?^*? 1 ' where Ta are the 

three generators of 0(3). In the representation where 

1 a0 — lt aal3, 

3 

This leads to the conserved (Noether) current: 

f i+& = ite ahc ^l^f. (13) 

(See Methods for its derivation.) One can then define a 
"conductivity" which is the linear response of this 0(3) 
current to an applied field. We shall show that this 0(3) 



conductivity has the classic marginal Fermi liquid be- 
haviour. 

Within the marginal Fermi liquid regime, the re- 
tains its marginal self-energy. The $ mode, while having 
no self-energy corrections at T = and uj = 0, acquires 
a Fermi liquid on-site self-energy because of the virtual 
fluctuations into the three (with the same diagram 
as Fig. 3, except that site i is the same as site j). In 
this regime then, our model has 2 distinct quasiparticle 
relaxation times: 

r 4 = -ImS^(ui) ocuj 2 + tt 2 T 2 (14) 
= — TmE^(cj) cx max[| lu |,T] (15) 

This suggests an intriguing link to the two relaxation 
time phenomenology observecLin the electrical and Hall 
conductivities of the cupratesE3, where electrical conduc- 
tivity is dominated by the slower relaxation rate, while 
the Hall conductivity is proportional to the product of the 
two relaxation times. Unfortunately, in our model, the 
conserved current j a does not include the $ fermion, thus 
transport quantities constructed from this 0(3) current 
do not reveal the slow, quadratic relaxation. But another 
conserved quantity in our model is the total energy. Fol- 
lowing the same strategy as for the 0(3) current, it can 
be shown that the conserved thermal current is just a 
sum of currents due to each of the i> and <j( a ). Now 
the propagators are diagonal in the $ or \P operators, 
thus the thermal conductivity proportional to the ther- 
mal current-current correlator will just be a sum of the 

relaxation times: k/T cx + j?-. Also, any mixed 
correlators (j a Q°) will be identically zero, where Q° is 
the thermal current due to In summary, it is unfortu- 
nately impossible to see the two relaxation times entering 
multiplicatively in any transport quantities of our model: 
the various conductivities derived from the 0(3) current 
will only depend on vfW, whereas the thermal conduc- 
tivity will be dominated by the largest relaxation time: 
that of the $ fermion. 

All of our previous considerations assume that the sys- 
tem does not develop into a broken symmetry state at 
low temperatures. In fact, as our Hamiltonian is a gen- 
eralisation of the Hubbard model at half-filling, it is per- 
haps not surprising that it displays similar magnetic or 
charge ordering instabilities, due to Fermi surface nest- 
ing. The main qualitative difference is the presence of a 
large marginal Fermi Liquid regime in the T — U phase 
diagram. 

The plan of the paper is as follows: in Section III, we 
set out the formalism of Dynamical Mean Field Theory 
for solving the lattice model in the weak coupling regime. 
Section IV presents the results and Section V discusses 
the lattice coherence scale, and the relationship to the 
two-channel Kondo lattice. We also discuss the low tem- 
perature phase of the lattice when the Fermi surface has 
strong nesting instability. Finally, we touch on the diffi- 
culties this model faces in modelling the cuprates. 
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III. METHOD 



We study the Hamiltonian: 



h = it j: $ 3+c $, + itj^ E 

j.c j,c a—1 



(16) 



where each fermion is a canonical Majorana fermion: 
{^\ *5 b) } = <U%, $,} = Sij, and {$f», $,} = 0. 

To gain some insight into the properties of the model 
in the weak coupling limit, we use dynamical mean field 
theory (DMFT)Ea, a method suited to systems where the 
dominant interaction is on-site and spatial fluctuations 
are unimportant, and when the on-site temporal fluctu- 
ations at all energy scales are to be taken into account. 
The d — > oo limit (where DMFT is exact), requires the 
usual scaling t a = t a j\Jd , with t a ,U ~ 0(d a ). For a 
lattice of localised bound states <F,-, to = 0. 

The crucial simplification in DMFT is that the self en- 
ergies are fc-independent as any intersite diagram (such 
as in Fig. 3) is at most of order \j \Jd relative to on-site 
diagrams. Now the diagram of Fig. 3 is precisely that 
which causes the $ fermion to propagate: omitting it 
means that to remains at zero, and $i stays localised, 
strictly in infinite spatial dimension. In finite d, we will 
need to incorporate its effect, since propagating <5>i fields 
will lead to the destruction of the marginal scattering 
mechanism. A rigorous 1/d expansion appears to be 
formidable. Nevertheless, the essential effects of finite 
dimensions may be included by introducing a finite value 
for to, and treat it as a fixed parameter of the model. 
Defining to = to/y/d, we estimate to by calculating the 
diagram (Fig. 3) in finite d at T = 0. In Appendix A, we 
show that t /t = c(U/t) 2 /d, with t a = t for a = 1,2,3, 
and c is a small constant. Finally, define a zeroth com- 
ponent = <Fj to get: 



H 



ta 



/O) ,T,( a ) 



a—0 v 



C 4 J J J J 



ta = 



3,c 
t 

to 



a 
a 



1,2,3 
0. 



(17) 



We shall use both $j and \Vj interchangea 

Following standard procedures of DMFTo, we map 
the lattice problem to an effective single-site problem 
with the effective action: 

S eff = f dTdT'{j2^ a) {T)g-\r-r')^\r') 



-U [ dT(*WyW*W¥V)(T), 
JO 



(18) 



where Q a (t) is the dynamical mean field at the single site 
which includes time dependent influence of the rest of 



the lattice, and is not the original lattice non-interacting 
local Green's function. It plays the role analogous to the 
Weiss mean field for conventional mean field theory, and 
has to be determined self-consistently: see below. 

The effective single-site problem dressed Green's func- 
tion G is related to Q via: 

E a (iw n ) = ^- 1 ^„)-G- 1 (iw„), a = 0,1, 2, 3, (19) 

where the self-energies E a are calculated from S e ff. To 
lowest order in U, they are given by the following dia- 
grams in Fig. 4. 





FIG. 4. On-site self-energies. 



The effective single site dressed Green's function must 
then be related back to the k— average of the original 
lattice dressed Green's functiop_G a (fc, iuj), via the mean 
field self-consistency equationslia: 



G a (iuJn) 



d d k 
(2tt)° 



G a {k,iuj n 



where since in d — > oo limit, self-energies have no 
fc-dependence, G a (k,iw n ) = 1/Uw n - e a {k) - £ (iw„)). 
Doing the integration then givestS 

Q^^iuin) = iuj n + it a sgn(uj n ), a = 0,1,2,3 (20) 

for the Lorentzian density of state (DOS): D a (e) — 
t a /(ir(e 2 J rt 2 a )) (corresponding to infinite range hopping); 
or alternatively: 

g-\ioj n ) =iLj n ~t 2 a G a (iuj n ), a = 0,1, 2, 3 (21) 



for the semi-circular DOS D a (e) = ^y/l - (e/2t a ) 2 
(corresponding to nearest-neighbour hopping on a Bethe 
lattice). Equations 17-19 and 20 or 21, and Fig. 4 together 
define our DMFT. The Lorentzian DOS is tractable an- 
alytically as the self-consistency equations are decoupled 
from G; however this also means that the effect of the 
lattice enters rather trivially as just a renormalisation of 
the bandwidth. To check these results we also use the 
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semi-circular DOS where the self-consistency equations 
are sobbed computationally using iterated perturbation 
theoryla. 

One quantity of particular interest in the context of 
the marginal Fermi liquid is the optical conductivity. As 
mentioned in Section II, the Hamiltonian of Eqn.17 does 
not conserve particle number unless to = t, thus the or- 
dinary electrical current proportional to the particle cur- 
rent is not useful. We can however generalise the concept 
of the optical conductivity to this model, using the fact 
that the total isospinB S a = "rM^T^^ is con- 
served, where T"^ = ie aa p are the 0(3) generators. Com- 
bining the continuity equation: id T S^+^2 & (j^ +a ~ j\ ) = 
(a are the unit lattice vectors), and the equation of mo- 
tion d T S\ = [H, Sf] leads to the conserved current: 



3i 



(22) 



This is the Noether current associated with the 0(3) 
symmetry of our Hamiltonian. We can then introduce 
a vector potential field A'°'(; 



( a ) A( a ) 



(Ay 1 a 



■A 



(ah 



in d 



dimensional space, coupled to the electrons as follows: 



tf^i^^exp jz J 



dl- A (a) T a 



(23) 



a0 



where the line integral goes from site i to i + c, and a 
summation is implied over dummy indices. Since there is 
isotropy for a = 1, 2, 3, we need only study the response 
to the the A^ component j x X (w) = J2 y a p±{u)A^ (w), 
which can be described by a Kubo formulaEd: 



-U(q, iuj n ) 



9=0 



U(q, iw n ) = - / dre^iTrjt (q, r) Jx (q, 0)) . 
Jo 

In the d — > oo limit, the_absence of vertex corrections to 
the conductivity bubblec_l permits us to write: 

o(v m ) = ^T^(v k J 2 G a (k,iuj n +iu m )G a (k,iw n ). 



As usual, at temperatures much lower than the band- 
width, doing the Matsubara sum leads to a function 
peaking largely near hp, and we replace J2k( v k m ) 2 by 
n/m J de. Doing the energy integral and analytically con- 
tinuing to real frequencies: 



a(y + id) 



did 



f(U-) - f(u+) 



x [ v -(Z R (Lo + )-X R (uj-))+i(T(Lu + )+r(uj-))} 1 (24) 
where E(w ± i0 + ) = =p iT(u>) and lo± = lo ± v/2. 



IV. RESULTS 

In the d — > oo limit (ie. to = 0), Go is the same as the 
single impurity model, by Eqn 19 or 20, so the bound- 
states described by are localised, with a self-energy 
that has a Fermi liquid form: 



£4,(0;+) 



(u/ + (*T)')}, (25) 



where uj + = oj + i0 + and Nq is the DOS at the Fermi 
surface. We give a brief derivation of this result in Ap- 
pendix B. For 5 , ,- a \ a = 1, 2, 3, the mean field Q a of the 
effective single site problem is of the same form as the 
single impurity model bare Green's function, using the 
Lorentzian DOS for the d — ► oo lattice. This is due to 
the DOS in the effective problem being smooth at Fermi 
energy, and for T -C t, equaling a constant No, just as 
in the single impurity modeO We obtain the marginal 
self-energy: 



£*(w+) = -{UNofu: 





f to \ 


7lT" 


In--* 






T 







(26) 



where A is a cut-off proportional to t, * is the Digamma 
function, 7 ~ 0.6 is the Euler constant (Appendix B). 
This has the following limiting behavior 



S*(w+) = [UNof x 



v M 
(-wln4ff 



iirT), 



(c^<T) 
(27) 



At finite d, the lattice coherence energy scale fo gen- 
erated from the diagram of Fig. 3 becomes finite, with 
to ~ U 2 /t. Marginal Fermi liquid behaviour will now per- 
sist so long as to < T <C t . At lower temperatures T <C 
t < t the three-body bound-states begin to propagate 
coherently, causing a cross-over to Fermi liquid behavior. 
This is borne out by analytical calculations; here we illus- 
trate with computational results (using the semi-circular 
DOS) in Fig. 5 showing the effective quasiparticle scat- 
tering rate T(u>) = LoRea(uj)/ Ima(u>) ~ — ImY< a (u) + ), 
where cr(u)) is the optical conductivity defined in Meth- 
ods. In Fig. 6 we plot the the resistivity p(T) showing 
the large linear T regime at weak coupling, and the inset 
shows the crossover to the T 2 Fermi Liquid regime. 



V. DISCUSSION AND CONCLUSION 

In this paper we have shown how a lattice of three- 
body bound-states induces marginal Fermi liquid behav- 
ior above a lattice coherence temperature t Q where a 
Fermi liquid forms. Since to ~ U 2 /t, a substantial win- 
dow in temperature for marginal Fermi Liquid behaviour 
exists only for small U. The emergence of this lattice co- 
herence energy is expected to be quite general: whenever 
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FIG. 6. Plot of resistivity vs.T. Inset shows the very low 
temperature crossover to T 2 behaviour: the y-axis is in the 
same units as the big plot. 



a localised mode is allowed to interact with conduction 
electrons, it will be difficult to prevent hybridizations be- 
tween fields at different sites. These effects however, are 
missed in a strict d — > oo calculation. 

This brings us onto the question of the relationship 
between our lattice model and the two-channel Kondo 
lattice model. The single-impurity version of the HamiL. 
tonian of Eqn.ll (with U < 0) was originally derivedEHl 
in the_cpntext of the single-impurity two-channel Kondo 
modeO. That derivation took advantage of spin-charge 
separationlij to throw away uncoupled (charge) degrees 
of freedom, and it has been shown via bosonization to 
be exactly equivalent to the original modelO. Unfortu- 
nately there is no such relation between the two-channel 
Kondo lattice and the Majorana lattice considered here. 

In the above calculations, we have assumed that the 
marginal Fermi Liquid state is unstable only to the Fermi 
Liquid state at low T. In fact, just as in the Hubbard 
model , a mean field calculation indicates that for a Fermi 
surface with a strong nesting instability (for example, 
nearest neighbour hopping in a hypercubic lattice), there 
is a phase transition to antifcrromagnetic order (for U > 
0). The order parameter is a vector that reflects the 0(3) 
symmetry of the model: 



— e iQ-Xj 



ie 



Vj/(0)vr,(a) _ _ ( I x i[. 



(28) 



where Q = (jr, . . . , 7r) is the nesting vector. From the di- 
vergence in the susceptibility, we find that at weak cou- 
pling, T c < to, except when to = 0: 



A 



exp 



In to A /(In i A) 2 1 + k/t 



2(UN y 



(29) 



where A is a cut-off (A < t). (Note that this reduces 
to the Hubbard model value when to = t.) At to = 0, 
T c is identical to the Hubbard case. For < to < t, T c 
is enhanced relative to the Hubbard case. Hence a re- 
gion of Fermi liquid phase separates the low temperature 
antifcrromagnetic phase from the marginal Fermi Liquid 
regime: Fig. 7. 

There are further similarities to the Hubbard model 
at half-filling. Both of the 50(4) and 0(3) models are 
invariant under U — > —U and * (0) — > The lat- 

ter map corresponds to a particle-hole transformation for 
the down spin only: C| <-> cf. This implies that in going 
from the positive U model to the negativCjXj model, mag- 
netic ordering turns into charge ordering^ . It can also 
be shown that neither model mixes charge and magnetic 
ordering. Further, both models reduce to the Heisen- 
berg antiferromagnet as U — > oo. Thus our model has 
very similar properties to the half-filled Hubbard model, 
except for the marginal Fermi Liquid phase at weak cou- 
pling. 

What insight docs our model bring towards the under- 
standing of the marginal Fermi liquid behaviour in the 
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FIG. 7. Schematic phase diagram of the Majorana lat- 
tice model. The dimensionless coupling is g = UNo where 
No <x l/t is the density of states at the Fermi surface. 
to — ctg 2 /d. MFL is the marginal Fermi liquid phase, FL 
is the Fermi liquid phase. AF is the antiferromagnetic phase, 
PI is the paramagnetic insulating phase. At weak coupling, T c 
goes as exp( — 1/g), while at strong coupling, it goes as t 2 /U. 
We have no computation for strong coupling for T > T c , thus 
we do not know the continuation of the T = to line. However, 
we might expect that when U is of order t, T c and to will also 
be of order t. One scenario is then the T = to line meets and 
ends at the T — T c line around to. Also, at strong coupling, 
the paramagnetic insulator crosses over to a metal at T ~ U 
(in analogy to the Hubbard model), and in the weak coupling 
regime, there will be a similar high temperature cross-over 
from marginal Fermi Liquid to an incoherent metallic phase: 
these cross-overs are indicated with a dotted line. 



cuprates? While our model does provide a simple lattice 
realisation of a marginal Fermi liquid, it unfortunately 
suffers from a number of defects: 

• It has a wide window of marginal Fermi liquid 
behaviour only for small coupling. The cuprates 
are believed to be in the strong coupling regimecj, 
but at strong coupling, our system has charge or 
magnetic instabilities. Related to this is the fact 
that in the cuprates, the inelastic scattering rate 
r = romaxfu;, T] has To/i a constant of order one, 
whereas our model has To/t proportional to the 
coupling squared. 

• The model needs to be at half filling: upon doping, 
a chemical potential term /i^ )^ 3 ) + vfrWvfrO 2 )) 
leads to a width A oc /x 2 for with Fermi liquid 
properties when T < A. This seems to require 
fine tuning, contradicting the rather robust linear- 
T resistivity observed even in underdoped systems 
(above the "spin gap" scale). 

• Despite the presence of two relaxation rates in the 
system, transport quantities will not involve a mul- 
tiplicative combination of the <f> and '5 relaxation 
rates, as is postulated in the two-relaxation-times 
phenomenology for the cuprates. (See Section II.) 

In conclusion we have demonstrated the persistence 
of marginal Fermi Liquid behaviour at weak coupling in 
a toy model of a marginal Fermi Liquid in an infinite 
dimensional lattice. For finite d, the lattice coherence 
energy cuts off marginal Fermi Liquid behaviour and the 
system reverts to a Fermi Liquid at low temperatures. 
Since this cut-off grows with the coupling, there will be 
no marginal Fermi Liquid regime at strong coupling. It 
remains to be seen if a strong coupling marginal Fermi 
Liquid exists in any finite dimensions. 
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VI. APPENDIX A: CALCULATION OF THE 
EFFECTIVE BANDWIDTH t 



We want to estimate in finite dimensions the effective 
kinetic energy ita J2i c ^i+c^f^ fr° m the zero-frequency 
part of the self-energy, as depicted in Fig. 3, to lowest 
order in the coupling: 



ito = U 2 



rfr[G 4 (r)] s 



(30) 
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where G$(t) is the bare propagator of ^ a \a = 1,2,3), 
for nearest neighbour sites, taken here to be in the x 
direction. In k, ui space, G(k,iu> n ) = [iu) n — eg] -1 - For 

simplicity, take = — 2f^^ =1 sin(fcj), appropriate for 
the hypercubic nearest neighbour dispersion. (We expect 
that as T — > 0, the exact shape for the band does not 
matter.) Let t = tj y/d and to = to/y/d, as required for a 
proper scaling of the kinetic energy term in large d limit. 
Then G x (t) will be of order 1/y/d, and as mentioned in 
the Methods section, to ~ 0(d~ 3 / 2 ), ie., 1/d down on 
the dispersion for the a = 1, 2, 3 components. Doing the 
standard Matsubara sum leads to: 

/d d k 
^d/(- e fc) ex P(« fc ' i - e fe T )- ( 31 ) 



In the zero temperature limit, the Fermi function become 
/( — e fc) ~~ ¥ ^( e fc» + where we have split up the disper- 
sion into the k x part and the other d — 1 part. Turning 
the d — 1 dimensional fc— integral into an energy integral 
gives: 



Gx(r) 



W)' 



,ik m -e k;c T 



deN(e)e- 



(32) 



where N(e) is the density of states in d — 1 dimensions. 

To make further progress, a flat density of states is 
used: N(e) = l/(4t) for |e| < 2t and zero otherwise. 
Thus, defining the dimensionless time variable s — 2rt, 



(33) 



where Jl(z) = || e -«sin(x)-Ma: j g B esse i f unc _ 

tion of the first kind of the first order. Thanks to the 
factor of 1/y/d inside its argument, the Bessel function 
asymptotically always grows more slowly than the decay 
due to the e~ s factor. (For |arg(z)| < tt and \z\ — > oo, 
Ji(z) — ► (^) 1 / 2 cos(z - 37r/4).) Thus, contributions to 
the integral in Eqn.30 are dominated by the regime when 
the Bessel function is at most of order one, allowing us to 
approximate: iJ\(ix) = —x/2 + 0(x 3 ) in Eqn.33, leading 
finally to: 



to 



El 

2t 



ds 



2s 2y/d 
U 2 



(34) 



2 7 -3 t d 3 / 2 ' 
with to of order d~ 3 l 2 as claimed. 



VII. APPENDIX B: CALCULATION OF THE 
MARGINAL SELF-ENERGY 



To order U 2 , the on-site self-energy for , a — 1, 2, 3 



is: 



s*(r) = u 2 g (r)g a (T)g a (T) 



(35) 



Go(t) = sgnkr)/2 is identical to the single impurity 
propagatorc3, since in the strict d — > oo limit, to = 0. 



For * 



(«) 



Qa(r) 



TV 

^ ioj n + it a sgn(u! n ) 



(36) 



where we have used Eqn.19 for Q a {iu) n ) for the Lorentzian 
DOS. As usual, turn the Matusbara sum into a contour 
integral and deform the contour onto the branch cut at 
the real axis to get: 



Q a {r) = -Im 



7T UJ -\ 



it a 



(37) 



(f(u>) is the Fermi function.) As we are interested in 
T <C t a , the integrand is dominated by small lj, and we 
approximate the denominator u> + it a « it a . Now the 
d— dimensional Lorentzian DOS at the Fermi surface is 
N = l/(7rt ), and thus GJa) is identical to that for the 
single impurity calculationE3: 



Qa(r) 



N irT 
sin(r7rT) 



(38) 



Note that this expression is accurate for -C t <C (3- 
Going to Matsubara frequencies u> n = (2n + l)nT: 



E*(iw n ) = V 



dr s 9 n i T ) 



N q ttT 
sin(r7rT)^ 

-i— (UN ) 2 I n (e) (39) 



where 



7r ~ e , sin(2n 
dx- 



1) 



(sin(x)) 2 



(40) 



We have put in a cut off e = nT/(2t a ) -C 1. Integrat- 
ing by parts twice, and using the tabulated integralEj: 
fg da;lnsin(a;) sin(2n + l)x = -^pr^^TT + m2 + 7 + 
*(n + 1/2)), we get: 



\2irTJ 



(41) 



ty(x) is the Digamma function, A = t a e 1-7 / (irT) and 7 ~ 
0.6 is the Euler constant. (We have expanded in e and 
kept only the terms up to e°.) Putting this all together 
and performing the analytic continuation iuj n — > u> + i0 + 
gives Eqn. 26. 

The on-site self-energy for $ is: 



^(T) = u 2 g a (T)g a (T)g a (T) 

Fourier transforming: 



(42) 
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£$(iu; n ) 
K n = 



iN (N UTrT) 2 K n , 
sin(2n + l)x 



(43) 



dx 



(sin(x)) 3 



Again, we have a cut-off e = TrT/(2t a ) <C 1. Integrat- 
ing by parts and expanding in e, K n = 2{2n + l)/e — 
2n(n + l)ir. This then leads to Eqn. 25 after analytic 
continuation. 
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